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Abstract: Diffraction limit is manifested in the loss of high spatial 
frequency information that results from decay of evanescent waves. As a 
result, conventional far-field optics yields no information about an object's 
subwavelength features. Here we propose a novel approach to recovering 
evanescent waves in the far field, thereby enabling subwavelength-resolved 
imaging and spatial spectroscopy. Our approach relies on shifting the 
frequency and the wave vector of near-field components via scattering on 
acoustic phonons. This process effectively removes the spatial frequency 
cut-off for unambiguous far field detection. This technique can be adapted 
for digital holography, making it possible to perform phase-sensitive 
subwavelength imaging. We discuss the implementation of such a system in 
the mid-IR and THz bands, with possible extension to other spectral regions. 
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1. Introduction 

Microscopic imaging is the oldest and one of the most important non-invasive analysis tech- 
niques. It has been immensely successful in uncovering the structure, composition, and dynam- 
ics of micro- and nanoscale chemical and biological samples. Contemporary investigations in 
the life sciences demand ever-increasing resolution in a variety of spectral bands, with much 
attention given to IR and THz. This task is made complicated by the diffraction limit, which 
sets a fundamental upper bound on the maximum spatial frequency conveyed by conventional 
refractive optics. 

Between the initial studies of optical resolution by Rayleigh and Abbe in the 19th century 
and the present day, a multitude of super-resolution systems and methods have been proposed 
and demonstrated (TJ. Extracting information beyond the diffraction limit remains an active 
area of research. One general strategy for achieving sub-diffraction-limited images originates 
in the fact that it is possible to increase the spatial bandwidth of an optical system by sacri- 
ficing certain other characteristics (e.g. field of view, temporal bandwidth, acquisition time). 
This idea was pioneered by Lukosz in his seminal 1967 paper \2\. Indeed, many contemporary 
super-resolution techniques can be viewed as implementations of the Lukosz approach Q]|2, 
including off-axis illumination [3|, structured illumination [4|, spatial frequency-shifting grat- 
ings [5] [6), and even near-field scanning microscopy QUI [9]. 

The idea of improving spatial resolution by using temporal degrees of freedom is partic- 
ularly appealing for two reasons. First, time multiplexing can help resolve ambiguities that 
arise when high spatial frequencies are scattered into the optical passband, as in the case of 
evanescent waves diffracting off a sub wavelength grating [6|. Second, in certain frequency 
bands outside the visible spectrum (e.g. far-IR or THz) it might be easier to manipulate sig- 
nals in the time domain than in the spatial frequency domain. The potential utility of time- 
domain (or temporal-frequency domain) for super-resolution has long been recognized; in 
fact, a frequency-multiplexing scheme involving conjugate moving gratings was described in 
Lukosz's original paper and demonstrated several decades later ITOlfTTI . However, owing 
to the complexity of the experimental setup, this scheme has not seen wide adoption. 
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Fig. 1 . Schematics of the proposed system. 



In the present paper, we propose a time-multiplexed super-resolution system that requires 
no moving parts and is based on coherent detection of a frequency-shifted signal. This scheme 
lends itself particularly well to super-resolved imaging in IR and THz. Our approach is based 
on a device that converts evanescent waves to propagating waves via diffraction on acoustic 
phonons. The scattered and frequency-shifted waves can be easily decoupled from the existing 
propagating spectrum that forms the regular diffraction-limited image (thus, the image is free 
from aliasing). With minimal processing, these shifted components can be used to distinguish 
subwavelength features. 

We will discuss two variations of this approach. Both rely on mixing the frequency-shifted 
fields scattered from the object with a reference wave, creating, at the detector, a beat note 
photocurrent which can be isolated through lock-in techniques. We will see that detection of 
high spatial frequencies is enabled by the spatial frequency offset of the scattered signal, and 
that a true super-resolved imaging configuration can be attained with an additional temporal 
frequency shift of the reference signal. 

2. Scattering from a phonon grating: general description 

The proposed super-resolution microscope/sensing system is shown in Fig. [T] The object is 
placed in the near field of an acousto-optic modulator (AOM) and illuminated with a plane 
wave from a mid-IR or THz source. Waves scattered from the object strike the phonon grating 
set up in the AOM by a running acoustic wave with frequency £2. Due to scattering on the 
phonons, the transverse wave vector k x of the incident radiation is shifted by integer multiples 
of the phonon wave vector q, while its corresponding frequency is shifted by integer multiples 
of £2. For a sufficiently large q, the evanescent components of the object's spatial spectrum 
(\k x \ > (o/c) can be scattered into the propagating waves with \kf x \ = \k x — q\ < ffl/c. The various 
spatial frequency components can be measured using a Fourier optics setup (e.g. a lens with a 
detector array in its focal plane). 

We model our system by considering a rectangular sound column (i.e. planar acoustic wave- 
fronts propagating in the x direction) interacting with a spectrum of incident plane waves in 



a dielectric medium. We neglect the diffraction of the sound field or multiple reflections, and 
assume weak interaction. Due to photoelastic effect, the sound field produces a sinusoidal mod- 
ulation of the dielectric permittivity fl2\ . 



e(x) = e + Aecos(gx — Qi), (1) 

which corresponds to a spatiotemporal volume grating. We may write the general form of the 
field inside the grating as a sum over the discrete diffracted orders, 

E = H A J ( z ) ex P[*'(** + jl) x ~ + Jtyt] • (2) 

j 

The scattered plane wave components A,-(z) are then governed by the Raman-Nath equa- 
tions Ifl3l . 



where k- 



Ae 

A'}(z)+e z Aj(z) = -^tf[A } -x{z) +A; + i(z)], (3) 

1/2 



, and C0j ~ CO. 

The amplitude of the j th diffracted order, Ay, is proportional to (Ae) 7 , with Ae/e <C 1, al- 
lowing to ignore higher order terms (j > 2). We can, furthermore, conclude that the amount of 
energy scattered into the shifted waves is small, thereby permitting to neglect the variation of 
th diffracted order Aq (the undepleted pump approximation) 1121 . We note that for propagating 
waves, this conclusion is valid insofar as there exists no Bragg matching between the incident 
and diffracted waves. Since the phonon wave vector q is a tunable parameter in our model, it 
is always possible to pick a range of q values to ensure minimal energy loss in the incident 
wave. For the evanescent waves, the undepleted pump approximation is justified by the small 
interaction length. 

Keeping terms up to first order in Eq. Q, we obtain: 

Ae 

Al(z)+k 2 z A ± {z) = -^co]A (z). (4) 

The scattering amplitudes of "upshifted" and "downshifted" waves can be obtained from this 
expression. Since the input field at spatial frequencies (k x =p q) contributes to the output field at 
spatial frequency k x , we may write: 

E ou t{kx) = [A- exp(/i2f ) +A + exp(— iGlt) +Ao] exp[/(^x - cot)], (5) 

with A± = t Efaikx^q), Aq = toEi n (k x ). Scattering coefficients t are, for the case of evanes- 
cent waves, given by 

± Ae (CO\ 2 1 Ae /on 2 1 

where k^r = y/e{co/c) 2 — (k x =fq) 2 , k — ik z — \/k\ — e(o/c) 2 , and fo is a Fresnel transmission 
coefficient. We see that the conversion of incident evanescent waves into propagating signals 
depends critically on the acoustooptic index contrast Ae and on the effective interaction length 
1 / K". We note, also, that Eq. |6]) describes also the generation of shifted spatial frequencies for 
the case where the incident wave is propagating (provided we make the association K = ik z ). In 
this case, the divergence of Eq. (|6| around q ss 0, as well as q ps ±2k x signifies the breakdown 
of the perturbative treatment of Eq. ([3]) due to the onset of Bragg-matching. 



From Eq. (|6]l we can estimate the diffraction efficiency of high spatial frequency input signal 
E- m (kf)as 

Ag (7) 
1 1 2kf n{\+n) 

with n = (the refractive index of the acoustic medium), and Ae °= \/F , the flux of acoustic 
energy per unit area. 

We thus obtain the amplitudes of the frequency-shifted waves. We define A± — t Ei n (k x ^fq), 
Aq = t()Ei n (k x ) (with linear coefficients to, t describing the generation of phonon-scattered 
and/or device transmission characteristics) and assume A,- 3> Aq ^> A , where A; is the detected 
amplitude of the illuminating wave Aje' k ° z . Averaging out the signal over the finite detector 
aperture and subtracting the background (which can be done electronically), we can write the 
intensity detected by the system of Fig.[T]as 

/out = |A | 2 + 2 (|A,A_ | 2 + |A,A + 1 2 ) 1/2 cos(i2f + y) . (8) 

The two terms in this equation can be decoupled using standard techniques: the DC term is 
isolated with the aid of a low-pass filter, while the term oscillating at the acoustic frequency 
Q. is recoverable using standard lock-in detection. For any given k x , this second term contains 
contributions from both A + = t + E- m {k x — q) and A_ = t~Ei a (k x + q). Although the coupling 
between these two quantities, together with the lack of phase information, makes it difficult to 
recover the spatial spectrum, the information collected can be used in detecting subwavelength 
morphological changes between different samples. 

3. Super-resolved fingerprinting: numerical simulation 

We now illustrate the ability of the proposed system to distinguished between subwavelength 
spatial features of different objects. In particular, we utilize Eq. ([8) to perform a comparison 
between the standard optical target (USAF test chart) and a modified target, where the label 
of every 6th line group has been randomly replaced. The first replacement corresponds to the 
last resolvable line group (A/2.5 line separation); the subsequent replacements correspond to 
halving the size of the line groups (A/5,. . ., A/40). We assume the measurement is performed by 
selecting an element of a photodetector array in the observation plane and using two orthogonal 
acoustic transducers to scan the acoustic wavevector within the range q x y £ [— 25 co/c , 25 co/c] . 
(We choose these values with the aim to improve by a factor of ~ 20 on the Abbe A /2 resolution 
limit, thereby collecting meaningful information about the A/40 line group.) 

In our computations, we assume the operating wavelength of 10 jim with germanium as the 
acoustic medium. We take Ae = 10~ 3 and restrict the magnitude of the acoustic wave vector q to 
25 co/c. Since for high spatial frequencies k™ « q, acoustic driving frequencies up to 8.75 GHz 
are required to retrieve k™ s» 25co/c. These parameters are within reach of modern ultrasonic 
transducers [ 14 1, as well as surface acoustic wave devices 1 15 1. 

It should be emphasized that any method that relies on digital processing of raw data can 
suffer form rapid - sometimes exponential fl6l - accumulation of noise. (Indeed, some pro- 
posed super-resolution methods even stipulate the need for an exponentially strong input signal 
to overcome noise and losses [ 17).) To show that this is not the case here, in our computations 
we add a normally-distributed random term to the AC amplitude of Eq. ([8]) in order to simulate 
noise in the system. Because SNR is expected to be lowest for maximum values of the acoustic 
wavevector q, we consider SNR=10 for q = 25 o/cQ Assuming a 20x20 element photodetec- 
tor array, we compute the signal given by Eq. ([8]) for the standard target, as well as the modified 
target [Fig. |2|a)]. Fig. |2jb) shows the result of subtracting the two datasets and performing an 

1 Acoustooptic diffraction efficiency, and hence the signal-to-noise ratio varies as \/q. 
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Fig. 2. (a) Optical test target and its modified version (inset). In the modified target, the 
"5" label of every column has been replaced by another digit, (b) Computed output of the 
system in the presence of noise (shown in grayscale) assuming a realistic, noisy detector 
with 400 active photocells. The modified optical target is superimposed for illustration 
purposes. The output of the system clearly identifies the location of every modified digit, 
even for regions far below the diffraction limit. 

inverse Fourier transform, with the resulting plot superimposed onto the modified optical target. 
Evidently, every change in the original image is manifested in this difference diagram. Further- 
more, it is largely localized in the vicinity of the actual changed pixels. It is possible to discern 
the difference signal even from the A/40 line group label. 

The ability to distinguish between fine spatial features of optical targets makes the system 
described above uniquely suited for identifying objects based on their subwavelength spatial 
features. As a result, it may find applications in fingerprinting and/or detection of chemical and 
biological structures. 

4. Super-resolved digital holography 

A straightforward modification of the setup described above not only allows to measure the 
"downshifted" A_ component directly, but also provides a method for retrieving phase infor- 
mation, making it possible to perform phase-contrast microscopy, and potentially enabling 3D 
imaging on subwavelength scales. 

To this end, a portion of the illuminating radiation is shifted in frequency by Q. b using 
a second AOM. Unlike the modulator that interacts with light scattered from the sample in 
the Raman-Nath regime 1121 . this second AOM utilizes an appropriately oriented and longer 
cell to produce Bragg scattering. This results in a strong optical signal at frequency CO + £l b , 
\Ab\exp[i(kb ■ r — (co + Q. b )t)], which is projected onto the detector [see Fig. [5Ja)] . Interfer- 
ence between the two optical signals produces beat note photocurrents with frequencies i2, Q. b , 

Iout(k x ) — \Eiexp(ikoz) +A b exp(ik ■ r) + [A~ exp(/£2f ) +A + exp(— iOt) +Aq] exp(ik • r)| 2 

= ...+2\A-A b \cos[{£l b +Q.)t +A<t>-]+2\A + A h \cos[{Q. b -Q.)t +A3> + ] +..., (9) 

where AO* = (k b — k) ■ r — <j) is the phase difference between the signal from the Bragg cell, 
\A b \ evp(ik b ■ r), and the Raman-Nath-scattered signal A ± exp(/A:- r) = |A ± |exp[/(0 ± +k-r)}. 




Fig. 3. (a) Schematics of the proposed system. Note the Bragg-shifted reference beam that 
aids in providing phase information, (b) The computed output of the system with optical 
test target as the object in the presence of noise. 



Of special interest is the component at frequency i2 + £lb, which carries the high spatial 
frequency information contained in its modulus and its phase A4>~ ~ (k x — k x )x — (j>~ . Both of 
these quantities can be retrieved using lock-in techniques. To produce the lock-in reference, the 
RF signals driving the two acoustic cells can be mixed using a nonlinear element (e.g. a diode) 
and appropriately filtered to produce the sum frequency. As a result, complete information can 
be obtained about the complex high spatial frequency Fourier component A~, from which it is 
straightforward to deduce the field Ei n (k x + q). By collecting data from multiple CCD pixels, 
as well as by varying the acoustic wave vector q, information can be collected about the entire 
spatial spectrum of the object. The data can then be digitally processed to produce a spatial- 
domain image containing subwavelength details, as well as phase contrast. 

Because the Bragg-shifted signal we use to decouple the A + and A_ terms serves as a refer- 
ence needed to record phase information, and because the image is reconstructed digitally, our 
technique belongs in the category of digital Fourier holography (DFH) |[T8l [191 120| . However, 
the proposed coherent detection method is different from traditional DFH setups in that it effec- 
tively converts sample spatial frequencies to temporal ones. In conventional holography, care 
has to be taken to isolate the target signal both in real and Fourier space. This translates into 
limitation on the field of view, as well as maximum attainable resolution ll20ll . The requirement 
that the CCD pixel spacing must allow for imaging the reference wave fringes further limits 
the resolution. As a result, simple digital holography setups suffer from very narrow spatial 
frequency passbands (N.A.^0. 1). Synthetic aperture techniques based on gratings 12111221 151 
have been successfully shown to increase the effective numerical aperture, as has heterodyne 
detection ll23l . Our method combines the benefits of these approaches: dynamic acoustic grat- 
ing allows to scan the spatial frequency space while always remaining in the CCD's passband, 
while the temporal frequency shifts multiplex the data, effectively increasing the number of 
information channels in a given spatial frequency band. 

We simulate the performance of the frequency-shifted digital holography system by first 
using Eq. |9]l to compute the response of the system to a calibration signal having unit amplitude 
for all spatial frequencies. In practice, such calibration signal might be generated by placing a 
point source in the vicinity of the AOM. Eq. (|9]l also provides the effective amplitude and phase 
transfer functions that allow to determine the detected signal for a given input field distribution. 
Gaussian noise is added to simulate spurious signals in the system. The input signal can then be 
obtained by dividing out the calibration quantities. In Fig. [3jb) we plot the simulated retrieved 



field magnitude. Zooming in on the central part of the test patten (figure inset) it is evident that 
every line group is distinctly resolved, suggesting that the effective numerical aperture is > 1 
(due to reconstruction of evanescent waves). 

5. Potential improvements of the proposed system and extension to higher frequencies 

As demonstrated earlier, the resolution of the proposed systems will depend on many factors, 
including the maximum attainable frequency shift, integration time, acoustooptic index con- 
trast, and success in minimizing detector noise, laser linewidth, and speckle. In addition to 
resolution improvement, there exist many possible ways to enhance the functionality and the 
performance of the proposed devices. For instance, because the phase information is preserved, 
the full complex field in the object plane can be reconstructed, potentially enabling 3D imag- 
ing via, e.g., phase-shifting interferometry 11241 . On the performance front, the sensitivity of 
the device may be improved with a subwavelength layer of highly doped semiconductor at 
the front AOM facet. When the dielectric constant of this layer is equal to —1, the evanescent 
fields are strongly enhanced due to resonant coupling to surface plasmons in the doped layer (a 
phenomenon known as "poor-man's superlensing" [25 1), leading to better SNR at the detector. 
Another possible way to improve the scattering efficiency of evanescent waves is placing the 
sample directly in the path of an acoustic wave, for instance, by running the wave through a 
microchannel containing objects to be studied. This approach may find many applications in 
novel integrated biological/chemical detection devices. 

Finally, we comment on the possibility of extending the proposed approach to frequencies 
other than the mid-IR and THz bands discussed here. While implementing the system for lower 
frequencies is essentially trivial, near-IR and optical frequencies pose a challenge. Acoustic 
phonon energies in practical devices do not approach the values necessary to produce a sub- 
stantial wave vector shift in these spectral bands. However, the required shift in spatial and 
temporal frequencies can in principle be attained by replacing the acoustooptic medium with a 
nanostructured periodically moving grating. In this scenario, the spatial frequency shift is de- 
termined by the periodicity of the grating while the temporal frequency shift is determined by 
the speed of grating oscillation. These parameters can be adjusted independently to optimize 
performance. 

6. Conclusion 

We have proposed a system that enables detection of sub-diffraction-limited spatial spectrum 
components in the far field by utilizing scattering from an acoustic grating. This process works 
whenever the spatial frequencies of the object are comparable in scale to the acoustic wave vec- 
tor. In its simplest implementation, the system could aid in "fingerprinting" of samples based on 
their subwavelength spatial features. With the use of an additional Bragg-shifted reference sig- 
nal, it is also possible to recover the phase of the original optical signal. The proposed approach 
has the potential to greatly enhance the specificity of mid-IR and THz spectroscopy. 
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